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We study the relative contribution of partonic sub-processes to D meson production and D meson- 
triggered inclusive di-hadrons to lowest order in perturbative QCD. While gluon fusion dominates 
the creation of large angle DD pairs, charm on light parton scattering determines the yield of 
single inclusive D mesons. The distinctly different non-perturbative fragmentation of c quarks into 
D mesons versus the fragmentation of quarks and gluons into light hadrons results in a strong 
transverse momentum dependence of anticharm content of the away-side charm-triggered jet. In 
p+A reactions, we calculate and resum the coherent nuclear-enhanced power corrections from the 
final-state partonic scattering in the medium. We find that single and double inclusive open charm 
production can be suppressed as much as the yield of neutral pions from dynamical high-twist 
shadowing. Effects of energy loss in p+A collisions are also investigated phenomenologically and 
may lead to significantly weaker transverse momentum dependence of the nuclear attenuation. 
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I. INTRODUCTION 

A useful probe of the dense nuclear matter created in 
collisions of heavy nuclei at the relativistic heavy ion col- 
lider (RHIC) is one that is sensitive to dynamical scales 
and can be both cleanly measured experimentally and re- 
liably calculated theoretically. Because of color confine- 
ment, only hard probes, i.e. those with large momentum 
transfers, can be reliably calculated in the perturbation 
theory of Quantum Chromodynamics (QCD) |lj,|2j. Ex- 
perimental measurements of inclusive particle suppres- 
sion 0, Q are now able to test jet quenching theory out 
to transverse momenta of px ~ 20 GeV [5| ■ On the other 
hand, typical dynamical scales of the nuclear matter pro- 
duced in relativistic heavy ion collisions are on the order 
of hundreds of MeV, which is both much smaller than the 
scale of a hard probe and non-perturbative. Therefore, 
an ideal probe should be not only "hard" but also sensi- 
tive to this "soft" physics. Open charm production has 
a potential to satisfy these criteria because of the two 
distinctive scales of the open charm meson: the charm 
quark mass, m c ~ 1.5 GeV, a relatively hard scale, and 
the binding energy, ~ Mo — m c ~ hundreds MeV, which 
may be relevant to the fragmentation and dissociation of 
D mesons. 

Little quantitative understanding of medium effects ex- 
ists for heavy-flavor production, due largely to the fact 
that the relevant data is quite sparse at present 0, 
and inferred from the modified cross sections for non- 
photonic electrons. The experimental situation is rapidly 
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changing, however. In addition to introducing new mass 
scales, the advantages of using heavy quarks as probes 
also arise from the fact that once formed, heavy mesons 
live much longer (~ 10~ n sec) than the duration of the 
QGP (~ 10~ 23 sec), travel macroscopic distances away 
from the creation point and can be measured directly. 
Such precision measurements are not only facilitated by 
the boost at forward rapidities at RHIC and the LHC but 
also impacted by the non-negligible cold nuclear matter 
effects in this kinematic region j^, 0, El El • 

With this in mind, we extend the study of Refs. |ToL 
IT^ | to hadronic collisions with heavy-quark final states 
and consider single- and double-inclusive open charm 
production. Because of the quantitative importance 
of multiple scattering effects for inferring properties of 
the nuclear matter, the models of in-medium interac- 
tions should be subjected to experimental verification for 
heavy-flavor production. One may begin to gain con- 
fidence by applying the theoretical description of such 
effects to heavy ion collisions in simpler situations, such 
as proton-nucleus (p+A) scattering. In the this case, 
plasma properties are not an issue and the main medium 
effect is, therefore, due to the multiple interactions in cold 
nuclei. The medium effects that we plan to investigate 
are also present in nucleus-nucleus (A+A) reactions dur- 
ing the interaction time Tint. — ^-Ra/i ^ T eq- ^ T QGP 
and cannot be neglected. Here Ra is the nuclear radius, 
7 is the Lorentz gamma factor of the nucleus, r eg . is the 
equilibration time and tqg p is the lifetime of the plasma. 

Multiple scattering affects the nuclear cross sections 
in various ways depending on how parton propagation 
is modified by the medium as the partons enter into, 
and emerge from the hard scattering. At one extreme, 
when the scattering from the medium is largely inco- 
herent, the parton modification is dominated by trans- 
verse momentum broadening 8]. At the other extreme, 
when the longitudinal momentum transfer is small com- 
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pared to the length of the path of the parton as it 
propagates through the nucleus, the scattering becomes 
coherent which can lead to attenuation, or shadowing. 
The coherent limit is described differently in different 
approaches 0, 0, ^j. In our work its effects are 
calculated in terms of nuclear-enhanced power correc- 
tions 0, 0, El . We have shown that these could ac- 
count for a large portion or all of nuclear shadowing ob- 
served at small Bjorkcn xb in lcpton- nucleus deep in- 
elastic scattering (DIS). They also contribute to the siz- 
able suppression of forward rapidity single- and double- 
inclusive hadron production at RHIC 0, 0|. As the 
heavy quark introduces a new mass scale, the depen- 
dence of the nuclear-size-enhanced power corrections on 
this scale needs to be elucidated. 

Inelastic interactions of partons propagating in nuclear 
matter lead to energy loss in both the incoherent 0] 
and the coherent [20l l2l| regimes. While jet quenching 
effects are best established in A+ A reactions HJQLOg, it 
was recently shown that energy loss can lead to hadron 
suppression in d+Au collisions at forward rapidities at 
RHIC 22] . In a preliminary study [2^ we approximated 
the effects of such energy loss by a rapidity shift. We 
will demonstrate that partonic implementation of energy 
loss, suggested in [gj], gives qualitatively similar results 
but also allows extraction of the fraction of the energy 
lost by an average parton in cold nuclear matter. Finally, 
we point out that the multi-parton dynamics discussed 
in this work will become even more important in p+A, 
and consequently A+A, collisions at the LHC. 

Our paper is organized as follows: In Section II we 
outline the theoretical framework for the calculation of 
observables in hadronic collisions in which elastic, inelas- 
tic and coherent multiple scattering effects in hot and 
cold nuclear matter can be naturally incorporated. In 
Section III we calculate to lowest order (LO) the con- 
tribution of the various partonic sub-processes to single 
inclusive D meson production and D meson-triggered in- 
clusive di-hadron yields in "elementary" nucleon-nucleon 
interactions. The light hadron/charm meson fraction in 
the away-side charm-triggered jet is shown to be sensi- 
tive to the transverse momentum of the associated par- 
ticle. The effect of nuclear-enhanced power corrections 
and their physical interpretation is discussed in Section 
IV. The suppression of both D meson production cross 
sections and of D meson-triggered correlations at forward 
rapidity at RHIC are also calculated. We discuss in Sec- 
tion V the effects of energy loss in cold nuclear matter and 
demonstrate that the resulting nuclear suppression has a 
different transverse momentum behavior when compared 
to dynamical high-twist shadowing. Our conclusions are 
given in Section VI. In Appendix A we explicitly calculate 
the LO parton scattering cross sections that contribute to 
single and double inclusive charm production. Appendix 
B contains the details of the factorization of the quark 
and gluon coherent scattering from the non-perturbative 
high-twist matrix elements. We also demonstrate how 
to reduce these high-twist matrix elements to the scale 




FIG. 1: Schematic representation of the initial- and final- 
state elastic, inelastic and coherent multiple parton scattering 
that give the nuclear ^"^-enhanced corrections similar to 
Eq. J2J for proton-nucleus collisions in the nuclear rest frame. 



of power corrections per individual nucleon enhanced by 
the nuclear size. 



II. THEORETICAL FRAMEWORK 

The interpretation of the plethora of results from past, 
present and future heavy ion experiments necessitates a 
reliable theoretical framework for the calculation of a va- 
riety of moderate and large transverse momentum ob- 
servables, ipip ^> Aq CD , that can be systematically ex- 
tended to incorporate corrections arising from the many 
body nuclear dynamics. The perturbative QCD factor- 
ization approach pj provides such a tool and expresses 
the observable hadronic cross sections as 

^hadronic — ^"partonic %ji /^rj f^f i: f^dj^) 

® jn^/fti^'woj ® In^ , : , /'... 'j • (i) 

where <g) denotes the standard convolution over the inter- 
nal kinematic variables in the reaction. In Eq. fifi 
and fj,dj are the renormalization, factorization and frag- 
mentation scales, respectively. (pi^^Xi, fif i) is the dis- 
tribution function (PDF) of parton in the hadron hi 
and D} lj /j(zj, /idj) is the fragmentation or decay func- 
tion (FF) of parton "j" into hadron hj. Factorization 
not only separates the short- and long-distance QCD dy- 
namics but implies universality of the PDFs and FFs and 
infrared safety of the hard scattering partonic cross sec- 
tions. 

In the nuclear environment, corrections in the basic 
perturbative formulas arise from the possible nuclear de- 
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Type of scattering Transverse momentum dependence of the nuclear effect 



Elastic (incoherent) 



Inelastic (final-state) 



Inelastic (initial- state) 



Coherent (t— channel) 
(u— channel) 



Coherent (s— channel) 



Inclusive 1-particle production: 
at high pt 

Inclusive 2-particle production: 

Inclusive 1-particle production: 
the kinematic bounds 
Inclusive 2-particle production: 
Pt correlations 

Inclusive 1-particle production: 
the kinematic bounds 
Inclusive 2-particle production: 

Inclusive 1-particle production: 
Pt dependence 

Inclusive 2-particle production: 

Inclusive 1-particle production: 
Pt dependence 

Inclusive 2-particle production: 



Cronin effect, enhancement at moderate pr, disappears 
di-jet acoplanarity, broadening of away-side correlations 
suppression at all pt, weak pt dependence, amplified near 
suppression of high pt correlations, enhancement of low 

suppression at all pt, weak pt dependence, amplified near 
same as inclusive 1 particle production 
suppression at low pt, disappears at high pt, pronounced 
same as inclusive 1 particle production 

enhancement at low pt, disappears at high px, pronounced 
same as inclusive 1 particle production 



TABLE I: Effect of elastic, inelastic and coherent multiple scattering on the transverse momentum dependence of single and 
double inclusive hadron production in the perturbative regime. 



pendence of PDFs and FFs as well as multiple soft parton 
scattering. These corrections may be enhanced by the 
nuclear size and the density of the quark-gluon plasma 
(QGP) in the deconfmed phase. In DIS, cold nuclear 
matter effects can be systematically organized in the 
framework of the factorization approach as follows: 

G hadron 

a™ (l + C^a s + C^al + cPal + ...)T^ 

(4) 

+ w ( 1 +^ 4) ^+^ 4) ^+4 4) ^+---) T(4) 

(6) 

+ ••• (2) 

In Eq. © T« are the twist "i" non-perturbative ma- 
trix elements. Q 2 is a typical virtuality, for example in 
hadron collisions we find Q 2 = —t, —u or — s. Standard 
higher order calculations at leading twist, such as NLO, 
NNLO correspond to the determination of the coefficient 
functions in Eq. J2J) . Soft in-medium interactions are 
controlled by new scales, such as the transverse momen- 
tum transfer squared per mean free path M 2 /A 0, Q and 
the scale of higher twist pernucleon£ 2 [HCDl. While we 
can write a formal expansion in a s for these quantities, 
it is usually subject to theoretical uncertainties and slow 



convergence of the perturbation series. For the purpose 
of this paper we consider the new scales to be properties 
of the medium and keep track only of the nuclear size 
A 1 / 3 , or L, enhancement. In DIS, high-twist shadow- 
ing corresponds to the identification and resummation of 
the nuclear-size-enhanced higher- twist terms, <rj T' 1 ' in 
Eq. @. We emphasize that in p+A collisions, illustrated 
in Fig. ^ neither initial- nor final-state interactions can 
be neglected since the coupling strength to the medium 
is qualitatively the same. 

From a phenomenological point of view, it is useful to 
classify the multiple parton interactions based on their 
experimentally observable effect on the hadronic spectra 
and the correlated yield of leading di-hadrons from away- 
side jets, or di-hadron correlations. Table [H] provides 
guidance, consistent with our current knowledge and re- 
sults from this work, to the correspondence between the 
nuclear modification to the cross sections measured in 
A+B collisions and the binary scaled p+p result, 



_ do-^g ' 1171 / dyi ■ ■ ■ dy n dp Tl ■ ■ ■ dp Tn 
AB ( n ab) da NN h " / d Vi--- dy n dp Tl ■ ■ ■ dp Tn 



In Eq. © (N A °g) can be evaluated from an optical 
Glauber model. 
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FIG. 2: Cross section for D° and D + charm meson produc- 
tion to LO in perturbative QCD at y/s = 1.96 TeV. Three ap- 
proximately equivalent pairs of a phenomenological i("-factor 
and the D meson fragmentation parameter r are shown. 



III. D MESON PRODUCTION AND 
CORRELATIONS IN LOWEST ORDER PQCD 

Nuclear-enhanced power corrections to the hadronic 
cross sections are sensitive to the partonic structure of 
the underlying hard scattering event. We first evaluate 
the differential cross sections for open charm production 
and open-charm-triggered di-hadron correlations in the 
Born approximation. Calculations to match the fixed or- 
der (FO) perturbative results to resummed logarithms 
of the type a™(afc h\pT/m c ) k up to next-to-leading log 
(NLL), n = 2, exist that treat flavor as "heavy" pij . 
based on the assumption 4> c /n{ x , [if) = 0. We include the 
charm contribution from the nucleon wave-function ex- 
plicitly since this approach leads to a faster convergence 
of the perturbation series [25| . Our goal is to investigate 
the fractional contribution of partonic sub-processes to 
the differential open charm cross section with power cor- 
rections in mind. We naturally limit our discussion to 
transverse momenta larger than the heavy quark mass 
PTt,PT 2 > m c . 



A. Partonic sub-processes for inclusive D meson 
production and D meson-triggered correlations 

The kinematics that underlies hadronic collisions is in- 
troduced by the following example of LO jet production 
in the collinear factorization approach: 



do 



jl (c)j2 (0 
NN 



dy c dydd 2 PT c d 2 PT d 



K 



NLO 



E 



S{Aip - tt)8(pt c -PT d 

PT d 



<t>a/N(Xa)(l>b/N(Xb) OC 2 — 2 
X "02 l M ab^cd\ 



(4) 



In Eq. Q S = (pa + Pb) 2 is the squared center of 
mass energy of the hadronic collision, x a = p^/pa + , 
x b = Pb Ipb~ are the lightcone momentum fractions of 
the incoming partons and 4> a /N{ x a), 4>b/N{ x b) are the dis- 
tribution function of partons a, b in the nucleon. In this 
work we use the CTEQ6.1 LO PDFs (26|. In the double 
collinear approximation the jets are produced back-to- 
back at LO in a s . The running of the strong coupling 
constant a s is also taken to lowest order. 

The power corrections to single- and double-inclusive 
hadron production are manifest only at small and moder- 
ate momentum transfers 0, 0] , where the effect of the 
heavy quark mass may not be negligible. PDFs, based 
on global QCD analysis, use massless evolution kernels to 
derive the Q 2 dependence of the heavy, c and 6, quark dis- 
tributions |23 . Systematic corrections of the type m 2 / Q 2 
have so far proven difficult to implement. To reconcile 
the use of parton distributions for charm quarks with 
m q = with a final state where m q ^ we calculate the 
relevant matrix elements in Appendix A. 

With rriTi = J p\. + mf , 
cone momentum fractions: 



c, d we identify the light - 



Xb 



1 
1 

7$ 



(m Tc e ya + m Td e Vd ) 



m T ,e 



(5) 
(6) 



From Eq. ijljl. after fragmentation into charm mesons, 
and integration over the away-side jets, we obtain the 
single inclusive D°, D + meson production cross section: 



aa NN 

dyid 2 p Tl 



= K 



NLO 



E 

abed 



dy 2 / dzx 
v Jv 



1 n , ,&a/N(Xa)<f>b/N(xb) a 2 — 



• (7) 



Here 



V 



{(si € [o,i])n(z 2 G [0,l])n(z x G [0,1])}; 

z\ = PTx/PTc and we have chosen the factorization, and 

renormalization scales /i/ = \i r = \jp\ + m 2 . Only 

charm quark fragmentation into D , D + mesons is con- 
sidered. 

In this paper we take D D / C (zi) for vector and pseudo- 
scalar mesons from 27]. The shape of the heavy quark 
fragmentation function is determined by the value of the 
non-perturbative parameter r. Smaller values of r cor- 
respond to harder fragmentation. The decay of the vec- 
tor states D *, D + * into D°, D + with the experimen- 
tally measured branching ratios is taken into account. 
To study the relation between the phenomenological K- 
factor in Eq. J7J) and the non-perturbative parameter 
r |27| in D meson fragmentation we compare the calcu- 
lated D° and D + cross sections to the Tevatron Run II 
= 1-96 TeV data 28] in Fig.0 We find that three dif- 
ferent combinations (-Kjvxo = 1<5, r = 0.1), (Knlo = 2, 



5 




4 6 
Pti [GeV] 

FIG. 3: The cross section for D° + D + charm meson pro- 
duction to LO in perturbative QCD at y/s = 200 GeV in 
p + p collisions at RHIC with Knlo = 2 — 3 and r = 0.2. 
The differential cross section is shown in four rapidity bins: 
y = 0, 1.25, 2.5 and 3.75. 
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FIG. 4: Fractional contribution of partonic sub-processes 
(z) to moderate and high pt x D° + D + meson production in 
y/s = 200 GeV p+p collisions at rapidities y = 0, 1.25, 2.5 and 
3.75. We considered (1) eg — > eg (solid line), (2) cq(q) — ► cq(q) 
(dashed line), (3) gg — > cc (dot-dashed line), (4) qq — > cc 
(double dot-dashed line) and (5) cc — > cc (dotted line). Note 
the dominance (except at the highest values of p:z\) of the 
eg — * eg channel. 



r = 0.2) and (Knlo — 3, r — 0.4) yield little difference in 
the single inclusive charm meson spectrum. We note that 
for the same choice of the parameter r — 0.1, a LO calcu- 
lation with standard 4> c /n(x, [if) ^ gives open charm 
cross sections similar to the ones from a NLO calculation 
that treats flavor as "heavy" . In this paper we adhere 
to the more conservative choice r = 0.2. Taking into ac- 
count an estimated ~ 50% uncertainty of the if -factors 
at the smaller y/s — 200 GeV, we show the differential 
D° + D + cross sections in p+p in Fig.Elfor several rapid- 
ity bins: y = 0, 1.25, 2.5 and 3.75. These are the base- 
line differential distributions relative to which the effects 
of power corrections and energy loss can be studied. The 
uncertainty of the choice at RHIC (K^lo — 2.5, r = 0.2) 
cancels in the nuclear modification ratios, Eq. J3J|. 

In hadronic collisions, heavy quark pair production 
proceeds via the "flavor creation" sub-processes, gluon- 
gluon fusion, gg — > cc, and light quark-antiquark anni- 
hilation, qq — > cc, at lowest order in the strong coupling 
constant a s . Since the gluon density is much larger than 
the quark and antiquark distributions at a small parton 
momentum fraction, x, the gluon fusion sub-process is 
believed to dominate charm production in hadronic col- 
lisions at collider energies. However, for inclusive single 
charm production, a single charm (or anticharm) quark 
can be created via the "flavor excitation" sub-processes, 
eg — * eg and cq — > cq. Although the charm parton dis- 
tribution is small, 4> c /n( x ) ^ 4>g/N( x )i such processes 
can be amplified by the matrix elements | M a ^ cc i\ 2 . In 
Eq. |0, over most of the available phase space 
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Li > 1 , 



t 




s 


u 




u 



I 2 »l, (8) 



where Li,L2 are numerically large. From the analysis in 
Appendix A, we find that the leading power behavior in 
Li of hard parton scattering for the first case in Eq. (JSJ 
is 



(\M g g^cc\ 2 ) ~ ~£i versus 



(9) 



<|M cg ^ C5 | 2 ) ~ ~L\ , (|M c ^ cg | 2 ) ~ 2L\ . (10) 

Clearly, there is large, cx L\, amplification of the scatter- 
ing rate in Eq. ifTUI) relative to Eq. © . 

The fractional weight of the partonic sub-processes to 
open charm production can be fully studied through the 
ratio: 



da 



dyid 2 p Tl I dyid 2 p Tl 



(11) 



Numerical results for the single inclusive charm cross sec- 
tions, Eq. |0, are shown in Fig. 01 at yfs = 200 GeV. 
We find that the dominant D meson production mecha- 
nism is the scattering of the heavy quark from the wave- 
function of the nucleon (or nucleus) on gluons, eg — * eg. 
Heavy on light quark scattering, cq(q) — > cq(q), also con- 
tributes significantly to the cross section. In these pro- 
cesses the charm anti-quarks from the nucleon wavefunc- 
tion end up in the fragmentation region near the beam 
rapidity. While gluon fusion, gg — > cc, is non-negligible 
it is also not the main channel for single inclusive open 
charm production. We have checked that i2 CT (pTi) i s 
practically independent of the interplay between Knlo 
and the hardness r of the c quark fragmentation into D 
mesons. 
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FIG. 5: Contribution of light hadrons and anticharm mesons 
to the fractional production cross section for pr t = 5 GeV 
D°+D+ meson-triggered away-side correlations for yi =y 2 = 
0. The insert illustrates the difference in the fragmentation 
of light quarks and gluons into hadrons versus that for the 
c quark into charm mesons. Three different values for the 
non-perturbative fragmentation parameter r = 0.1, 0.2 and 
0.4 were used. 



B. Hadron composition of D meson-triggered 
away-side jets 

We now proceed to charm meson-triggered di-hadron 
correlations in away-side jets, which have not been stud- 
ied previously. For pairs of hadrons with large open- 
ing angle between them, Aip — ipi — ip 2 ~ tt, includ- 
ing charm- anticharm and charm-hadron correlations, we 
find ESEl: 



dy 



= K NLO Y^ [ ^±D Dl/c { Zl )D h2/d {z 2 ) 
idy 2 dp Tl dpT 2 ^ d Jv Z\ 

(12) 



0a/N(Xa)<t>b/N(Xb) O 2 — 2 

X 7^ \Mab^cd\ 

Xa%b & 



where 



v = {{xi e [o, i]) n (x 2 e [0, i]) n { Zl e [o, 1] n (z 2 e [o, 1])} 

The fragmentation scale for light hadrons in Eq. 112|) 
is set to /x = pT d /z 2 and z 2 = z\pt 2 Ipti_- We use 
FFs from |3(j for , K ± and p(p) . In the calcula- 
tions that follow, we also include charm quark fragmen- 
tation into light hadrons. For the charm anti-quark, 
D D/c( z ) = D D/c (z). 

The expectation for non-trivial pt 2 dependence of 
charm-triggered away-side correlations is based on 
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FIG. 6: The double inclusive hadron production cross section 
at RHIC for a D° + D + trigger at p Tl = 5 GeV. Three differ- 
ent combinations of r and Kmlo are studied. Two different 
rapidity gaps y\ — yi = 0, 2 are shown for comparison. The 
characteristic "hump" represents hard c fragmentation into D 
mesons. 



the very different behavior with respect to z = 
phadron j ^parton Q f tne frag mentation functions of par- 
tons into light hadrons [3(| compared to those of heavy 
quarks into charm and beauty mesons |27|. The insert 
in Fig. [3] shows that while light hadrons favor soft de- 
cays of their parent partons, heavy quark fragmentation 
is very hard. For z ~ 0.6 — 0.9 there is more than an or- 
der of magnitude enhancement in the decay probabilities 
D D o/c(z) + D D+/c (z) relative to (l/N) Ei=i,Ar(Ar±/i + 
Dk± 1% + Dp{p)/i)i where i runs over the light and charm 
quarks, antiquarks and the gluon. 

Triggering on a D meson fixes the momentum of the 
charm quark much more reliably than does triggering 
on a light hadron [Tfl |2^| in studying away-side corre- 
lations. The non-perturbative fragmentation therefore 
controls the yields of di-hadrons versus the associated 
momentum px 2 and the abundances of the different par- 
ticle species. Fig.|S]shows our prediction for the hadronic 
composition of the D° + Z) + -triggered away-side jet. At 
transverse momenta significantly smaller than the trig- 
ger transverse momentum, pt 2 <§C pt ± , the away-side jet 
is dominated by pions, kaons and protons. At transverse 
momenta pr 1 ~ pr 2 , the away-side jet is expected to 
be dominated almost completely by D° and D~ mesons. 
Since there is little sensitivity to the choice of r, our 
prediction is robust and can be used as a test of the pro- 
duction mechanism of heavy quarks. 

Fig. shows the behavior of the Z?-meson-triggered 
di-hadron cross section, Eq. (|12|) . The differential spec- 
trum cannot be described by a single power law for 
Pt 2 — PTi ■ In the case of hard fragmentation, for ex- 
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ample r = 0.1, the double inclusive distribution can even 
be non-monotonic as a function of the associated particle 
momentum px 2 ■ The characteristic "hump" is associated 
with the c fragmentation into D° and D~ . This signa- 
ture is unique for heavy quark triggers and is seen to be 
independent of rapidity. For pion or unidentified hadron 
triggers one expects a smooth power law behavior as a 
function of pr 2 0,H3- We conclude, through the com- 
parison of Figs.0and[U that charm-triggered correlations 
provide a much more powerful handle on the hardness of 
D D / C (z) than does the single inclusive spectrum. 

In summary, even in p+p collisions at RHIC and the 
LHC, triggering on charm or beauty mesons can provide 
constraints on the heavy quark production and fragmen- 
tation mechanisms that are not accessible via the Teva- 
tron single inclusive data. Such measurements may al- 
ready be possible at RHIC using electron or muon trig- 
gers as proxies for heavy quarks |3lL l32j| . 



IV. NUCLEAR-ENHANCED POWER 
CORRECTIONS TO OPEN CHARM 
PRODUCTION 



In this section we focus on the effects of coherence 
that arise from the multiple soft scattering of the par- 
tons that participate in the hard interaction. The pro- 
cesses where such effects can be cleanly studied are deeply 
inelastic scattering (DIS), illustrated in Fig. 0a), and 
Drell-Yan. If the longitudinal momentum transfer along 
a fixed direction becomes small, the probe will inter- 
act simultaneously with more than one nucleon, due to 
the uncertainty principle. In a frame where the longi- 
tudinal size of the exchanged virtual meson is given by 
l c = 1/xP the Lorentz contracted nucleon has longitu- 
dinal size 2ro/7 = 2ro/(P/mjv). The critical value of 
the momentum fraction for the onset of coherence can be 
easily estimated to be xn — l/(2r mjv) ~ 0.1 [l2l Il3j. 
For x < 0.01, the interaction with all of the nucleons at 
the same impact parameter b is coherent, given the small 
transverse size A± = l/Q 2 , Q 2 > m 2 N , determined by the 
virtual probe. The multiple final-state interactions with 
these nucleons are mediated via the field of soft gluons 
and are formally suppressed by a power of the squared 
momentum transfer l/Q 2 for each additional coherent 
scattering. The resummed contributions form all powers 
in l/Q 2 have the effect of a dynamical generation of a 
quark or gluon mass as we shall see below. 

In the case of hadronic collisions there is indeed a sim- 
ilarity with the DIS dynamics in the final-state rescat- 
tering of the struck, small Xb parton from the nucleus, as 
shown in Fig. 0b) . Such nuclear-enhanced power correc- 
tions lead to suppression of the single and double inclu- 
sive hadron production cross sections at forward rapidity 
as long as the coherence criterion is satisfied [Kj . 



A. Physical interpretation 

The physical interpretation of the resummed nuclear- 
enhanced power corrections is most transparent in the ex- 
ample of neutrino-nucleus scattering because of the pres- 
ence of a physical mass scale M in the final state in some 
of the sub-processes 01 . Even in the absence of final- 
state interactions, the charged current s-quark to c-quark 
transition, s ± W ± — > c, requires rescaling of the value 
of Bjorken xb = Q 2 /2ra^v in the parton distribution 
functions 



xb 



x b 1 



AP 
7W 



(13) 



Since x represents the lightcone momentum fraction car- 
ried by the struck quark, Eq. l(T3|) expresses the con- 
version of the larger initial parton energy into the con- 
stituent mass of the heavy quark. In Ref. 13], using 
Dirac equations of motion, we showed that the physical 
mass and the dynamical power corrections commute and 
therefore the two effects add. The corresponding change 
in the value of Bjorken x is 



x = xb 



x = xb ( 1 



x B 1 



ejA 1 / 3 - 1) ap 

Q 2 + ~Q 2 



2 

711 dyn 


M 2 


Q 2 


+ Q 2 


%) 





L dyn 



(14) 



Our resummed higher-twist corrections rescale the mo- 
mentum fraction, xb, in Eq. I|14|) in a manner identical 
to the modification arising from the constituent heavy 
quark mass, Eq. ljT3)l. We identify m 2 dyn = ^{A 1 / 3 - 1) 
with the dynamical mass 

Elll 113 generated by prop- 
agation of the struck parton through the nuclear back- 
ground chromo-magnetic field. We remark that this ef- 
fect, which we calculated in QCD, is analogous to the 
generation of dynamical mass for electrons propagating 
in a strong electro-magnetic field (34|. For xb < 0.1, 





FIG. 7: Coherent multiple scattering of the struck parton 
in deeply inelastic scattering (a) and in proton-nucleus col- 
lisions (b) . Notice the difference between p + A and DIS in 
terms of the multiple initial- and final-state scattering of the 
incoming and outgoing partons ("a" and "c"). 
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Effective quark x 

Q = 2 GeV 

CTEQ6 LO parton distributions 



FIG. 8: Effects of rescaling of the momentum fraction x 
illustrated on the example of the lowest order CTEQ6 parton 
distribution functions for valence and sea u-quarks. Gluons 
rescatter in the final-state to lowest order only in hadronic 
collisions and are included for completeness. The effective 
parton distributions or partonic flux are reduced after the 
rescaling of the momentum fraction. We used Q 2 — 4 GeV 2 
and the large rescaling in x is for purposes of illustration only. 



due to coherence, 



m 



dyn 



has to be incorporated in the 



underlying kinematics of the hard scattering, Eq. (|14fl . 
For values of xb > 0.1 rescattering of the struck parton 
still occurs but the vector meson exchange is localized to 
one nucleon. The elastic final-state interactions vanish 
by unitarity upon integrating over the full phase space 
for the hadronic fragments, X. 

The mechanism behind the reduction of the deeply 
inelastic inclusive scattering cross section and the dif- 
ferential hadron production cross sections is illustrated 
schematically in Fig. [S] PDFs arc falling functions of the 
momentum fraction x: 



V (X) OC x- l ' 2+S v , 

<p a (x) oc x- 1+s ° , 

<j> s ( X ) CK X~ X+S ' , 



valence quarks , 

sea quarks , (15) 

gluons , 



where |Au, a ff | <C 1- The evaluation of the parton distri- 
bution function at an effectively larger value of x leads 
to a reduction of the cross section, — Act From 
Eq. 1151) it follows that the same dynamically generated 
parton mass for quarks apparently produces a different 
reduction in the sea and valence contributions to the 
DIS cross sections, — Arj soa ~ — 2Acr va i oncc . If a gluon 
scatters coherently in the final state, the generated dy- 
namical mass is Ca/Cf =9/4 larger than the dynami- 
cal mass for a quark due to the larger average squared 
color charge. Consequently, the observable "shadowing" 
of gluons is substantially larger than that for sea quarks, 
despite the very similar low-x behavior of the PDFs in 
Eq. (|T3|l . — A(Tgi uon ~ ~9/4Acr sca . We can, therefore, 



derive a hierarchy of dynamical high- twist shadowing, 

cr(x,Q 2 ) 

S v (x, Q 2 ) < S s (x, Q 2 ) < S g (x, Q 2 ) , (16) 

for cross sections with multiple scattering dominated by 
valence quarks, sea quarks or gluons, respectively. 



B. Power corrections in p+A reactions with heavy 
quark final states 

A detailed derivation of nuclear-enhanced power cor- 
rections in proton(deuteron)-nucleus collisions is given in 
Appendix B. In this subsection we summarize the main 
results. Unlike structure functions in DIS, hadronic col- 
lisions involve a more complicated momentum flow even 
in LO tree level diagrams. This requires isolating the Xb 
dependence of the integrand in Eqs. Q) and (|T2)l in the 
function 



Fab^cd(xb) 



(t>b/N(xb)\M ab ^cd\ 
Xb 



(17) 



Since we use the standard on-shell hard scattering ma- 
trix elements \M a b^ c d\ 2 reproduced in Appendix A we 
obtain an upper limit on the dynamical high-twist shad- 
owing. The scale of high-twist corrections on an indi- 
vidual nucleon £ 2 = 0.09 — 0.12 GeV 2 is fixed by the 
world's data on deep inelastic scattering on a wide range 
of nuclear targets [E|, including the shadowing in the 
structure functions, F%, and the enhancement in the lon- 
gitudinal structure function, Fl, relative to the result at 
leading twist. This is the same scale that was used to pre- 
dict the Q 2 -dependent suppression of the small x struc- 
ture functions in neutrino scattering on iron (Fe) targets 
measured by the NuTeV collaboration [3^. Since £ 2 is 
constrained from final-state quark scattering, we account 
for the representation of parton d in hadronic reactions 
by the color factor Cd on a process-by-process basis. The 
values are C q ^ = 1 and C g = Ca/Cf = 9/4 for a quark 
(antiquark) and a gluon, respectively. 

The forward rapidity particle production in p(d)+A 
collisions is analogous to DIS, and resumming the co- 
herent scattering with multiple nucleons is equivalent to 
a shift of the momentum fraction of the active parton 
from the nucleus which leads to a net suppression of the 
cross sections. The hard scattering scale is given here by 
i = q 2 = (x a pA — p c /zi) 2 , and for Eqs. Q and ljT2)> we 
derive the following nuclear modification: 



F a b~>cd(Xb) 

Fab — *cd 1 Xb 



1/3 



1) 



-t + m 2 



(18) 



We denote by A the atomic mass number of the nucleus 
even though it is labeled as hadron "£?" or "6" . The 
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FIG. 9: Multiple coherent scattering of the incoming and out- 
going partons in p+A reactions in the t-, s- and ^-channels, 
respectively. 



i-dependence of this shift in Eq. i|18|) indicates that the 
attenuation increases in the forward rapidity region but 
disappears at large transverse momenta. We note that 
the effect of the mass of the final-state parton, for exam- 
ple a charm quark, in Eq. (|18|l is precisely the same as in 
the DIS case pjj, Eq. ijg. 

While we have focused in detail on the power correc- 
tions in the ^-channel, it should be noted that the co- 
herent initial- and final-state interactions of the incom- 
ing parton from the proton, see the bottom diagrams 
in Fig. |(?JJ, can also be calculated. The results that we 
present here are analogous to Eq. 118fl . 



F a b^cd{Xb) 

Fab — >cd ( h 
F a b^cd{Xb) 

Fab — >cd I %b 



£2 (A l/3 _ y 



i + c a 



|VV3-1) 

+ ° c ~, 2 



(19) 



(20) 



In Eqs. (|19|l and (|20|) C a and C c are the numerical coeffi- 
cients analogous to Cd corresponding to partons a and 
c, respectively. By direct comparison of the nuclear- 
A 1 / 3 -enhanced power corrections in Eqs. 118(1. I|19|) an d 
(|20|l it is easy to see that in the forward rapidity region 
\t\ <C \s\,\u\ only DIS-like dynamical shadowing correc- 
tions are important. 

On the other hand, in processes with no strong final- 
state interactions, such as the Drell-Yan, coherent mul- 
tiple initial-state scattering can give only enhancement 
of the cross section, as shown in Eq. I(19J) . because of the 
rescaling of Xb to smaller versus larger values. Such result 
is also easily understandable since an incoming parton 
with large dynamically generated mass prior to the hard 
scatter has to carry a smaller fraction of the nucleon mo- 
mentum in order to produce particles in the final state 



with fixed momentum. Thus, initial-state multiple soft 
interactions without energy loss always lead to a Cronin 
effect 8] independently of whether they are treated as 
coherent or not. 

The derivations given here are critical to elucidating 
the dynamical origin of nuclear effects in high energy 
hadronic reactions. We have shown through explicit cal- 
culation that such effects are process dependent and may 
change their sign. Therefore, these are not factorizable 
as a part of the PDFs and FFs. 

The cross section for single and double inclusive hadron 
production can be calculated from Eqs. Q) and (|12H via 
the substitutions Eqs. p8|) . (|19fl and l|2(J|l in minimum 
bias reactions. We recall that the A 1 / 3 -cnhanced scale 
of higher twist, £ 2 = 0.12 GeV 2 , has the minimum bias 
geometry included. For specific centrality classes we ac- 
count for the varying nuclear thickness along the line of 
parton propagation as follows: 



A l/3 ^ A l/3_ 



T A (b) 



^^(fomiii.bias) 

where T A (b) = f™ p A (z, b)dz. 



(21) 



C. Numerical results 

To illustrate the similarities and differences between 
massless and massive final-state partons, we carry out 
a comparative study of the effect of power corrections 
on single and double inclusive ir° and D meson produc- 
tion. The left panels of Fig. show the suppression 
of the low and moderate px t neutral pion cross section 
relative to the binary scaled p+p result, Eq. ©, from 
high-twist shadowing. The nuclear modification factor is 
shown for two different forward rapidities, y\ = 1.25, 2.5, 
in *fs = 200 GeV d+Au collisions at RHIC. At trans- 
verse momenta pr ± = 1-5 GeV the suppression can be 
as large as 25% but disappears toward higher pt x due 
to the power law nature of the effect. RdAu{PT x ) also 
shows non-zero rapidity dependence with the effect in- 
creasing at forward rapidities. For comparison, we show 
the PHENIX measurement of the nuclear modification 
of muons coming from the decay of charged hadrons in 
mid-central d+Au reactions at yi = 1.8 is shown [l7| . 
Coherent power corrections cannot fully account for the 
nuclear suppression measured by experiment. The dis- 
crepancy becomes larger at higher transverse momenta. 
The bottom left panel of Fig. ^| shows the centrality 
dependence of coherent power corrections for light pi- 
ons. From Eq. I|21|l. we see that in central collisions the 
larger nuclear thickness generates correspondingly larger 
dynamical parton mass leading to larger suppression of 
the cross sections. 

The right panels of Fig. ^| show the calculated sup- 
pression of D° + D + mesons (and equivalently, D° + D~ 
mesons), in deuteron-gold collisions at RHIC. It should 
be noted that the nuclear modification is very similar 



1.2 - 

CD 

in 

ii 0.8 h 



0.6 - 



0.4 - 



0.2 





1.2 - 



I 

j? 0.8 

< 

n 



0.6 - 
0.4 - 

0.2 



"■"n — 1 — r -1- ! — 1 — r 

Minimum bias" 

Light pions (ji 




PHENIX 20-40%/60-88% 
at y = 1 .8 

J I I ■ I I L_L 



Centrality dependent 

y, = 2.5 




No E-loss 



"I — ' — T -1 — I — 1 I 1 

Minimum bias " 



D mesons (D +D + ) 



y, = 1.25 

. y, = 2-5 

• PHENIX prelim, p.", y = 1.6 
_i i i i i i i i i_ 



Centrality dependent 
y, = 2.5 




— - b = 3 fm 

b = 5.6 fm 

■ - ■ b = 6.8 fm 

i I . I . 



12 3 4 
Pti [GeV] 



12 3 4 
Pn [GeV] 



1 i 1 i — 1 — r -1- 1 — r 

1.2 p y 2 = Minimum bias — p 

o o 

Jt - 71 



S 0.6 - 



0.4 



1.2 - 



< 0.8 - 



0.6 - 



0.4 




No E-loss 

i ■ i L_ 



L_ 



I 



Centrality dependent 
y, = 2.5, y 2 = 



1 I 1 — I — 1 I 1 I — r 

^2 = Minimum bias 

(D°+D + ) - (D°+D 




— y,=2.5 _| 

— ^ = 1-25 

■ I L 



I 



Centrality dependent 
y, = 2.5, y 2 = 



-- b = 3fm 
— b = 5.6fm' 
■- b = 6.8fm 
J .l.l l_ 



1 2 3 4 5 
Pn = P T 2 t GeV l 



1 2 3 4 5 

P-T1 =P T 2[ GeV l 



FIG. 10: Top: suppression of the single inclusive hadron 
production rates in d+Au reactions versus pr for rapidities 
yi — 1.25 (smaller attenuation) and j/i = 2.5 (larger atten- 
uation). Bottom: impact parameter dependence of the cal- 
culated nuclear modification for central, 6 = 3 fm, minimum 
bias, fcmin.bias = 5.6 fm and peripheral, b = 6.8 fm, collisions. 
The scale of higher twist per nucleon £ 2 = 0.12 GeV 2 . Data is 
from PHENIX [TtI |36|1 . Theory does not include energy loss. 



FIG. 11: Top: suppression of double inclusive hadron pro- 
duction rates in d+Au reactions versus pr for trigger and 
associated hadrons of the same transverse momentum but dif- 
ferent rapidities: y\ = 1.25 (smaller attenuation) and yi — 2.5 
(larger attenuation) for the trigger and i/2 = for the associ- 
ated hadron. Bottom: impact parameter dependence of the 
calculated nuclear modification for central, 6 = 3 fm, mini- 
mum bias, fomin.bias = 5.6 fm and peripheral, b — 6.8 fm, colli- 
sions. The scale of higher twist per nucleon £ 2 = 0.12 GeV 2 . 
Theory does not include energy loss. 



to or slightly larger than that for ir°. The reason for 
this similarity is that in both cases the dominant channel 
of hadron production is via quarks scattering on gluons, 
which have the same (large) color singlet coupling to the 
medium, Cd — 9/4 in Eq. (JTSJ. In addition, the typical 
momentum fraction, z\, which enters in the determina- 
tion of the hard scale, t cx l/zi, is slightly larger for the D 
mesons. Preliminary PHENIX data on the modification 
of ^i~'s coming from the decay of heavy flavor at y = 1.6 
in cold nuclear matter is provided as a reference [36( . Al- 
though the error bars are large, D mesons seem to be 
suppressed as much as light hadrons. This attenuation 
effect is larger than predicted by power corrections alone. 

Our results for the dynamical shadowing effects in in- 
clusive two-particle measurements at large angle ~ 7r 
are given in Fig. 1111 We have chosen px 1 — Pt 2 anc l the 
same rapidity of the trigger hadron, yi, as in Fig.^| The 
associated hadron was chosen at rapidity yi — to facil- 
itate the study of di-jets with a large rapidity gap. The 
left panels of Fig. 1111 show that 7r° — ir° modification in 
the nuclear medium follows qualitatively the suppression 
of single inclusive pions as shown in Fig. 111)1 The right 
panels of Fig. ^2 illustrate the effects of coherent power 
corrections on (D° + D + ) — (D Q + D~) pair production. 
In this case we have a heavy c quark scattering coher- 
ently in the final state. At pr 1 = pt 2 — 1-5 GeV ~ m c , 



the effect of the heavy quark mass on the power correc- 
tions that we derived in Appendix B is noticeable, see 
Eq. Q1H[). At higher transverse momenta the mass effect 
is not relevant. Naively, we would expect that DD pair 
production would be less suppressed due to the quark 
versus the gluon color singlet coupling to the medium, 
C c i = 1 versus 9/4. However, in light hadron pair pro- 
duction the momentum fractions z\ , Z2 are significantly 
smaller |38j| than in single inclusive measurements and 
D meson production. Therefore the relevant hard scale, 
t, in DD pair production is smaller which amplifies the 
effect of the dynamical heavy quark mass. 



D. Discussion 

The nuclear-size-enhanced power corrections that en- 
ter this work are the result of a tree-level QCD calculation 
and represent the effect of multiple coherent scatterings 
of the outgoing partons with the gluon field of the nu- 
cleus. The strength of this interaction is specified by a 
universal parameter, £ 2 , which is related to the limit- 
ing soft gluon distribution of a nucleon in the nucleus. 
Once £ 2 is specified, the effect of coherent scattering for 



all partons, heavy and light quarks as well as gluons, is 
completely determined. In particular, when the coher- 
ence length, defined as the inverse of the longitudinal 
momentum transfer, becomes comparable to the nuclear 
size, these power corrections acting on partons in the fi- 
nal state suppress the hadron cross sections. 

In Ref. [lg, the effects of coherence are determined 
by multiple scattering, and this approach is therefore, 
in principle, quite similar to ours. The following points 
of difference should, however, be noted. Using a Green- 
function description of quark-gluon Fock states in com- 
bination with the Glauber-Gribov formulation of nuclear 
shadowing, it is shown there how multiple diffractive 
scattering of gluons suppresses the gluon distribution in 
a nucleus when the coherence length in gluon-dominated 
processes begins to exceed the inter-particle spacing, thus 
providing an estimate of gluon shadowing. The scale 
of the suppression is fixed by, among other things, the 
triple-pomeron coupling G3P — 3 mb/GeV 2 . The small 
triple pomeron coupling suggests a rather confined gluon 
cloud of valence quarks, and this weakens gluon shadow- 
ing and delays its onset. This gluon shadowing is also 
"leading twist", i.e., weakly dependent on Q 2 , vanishing 
logarithmically as Q 2 — > 00. 



V. ENERGY LOSS IN COLD NUCLEAR 
MATTER 

In Fig. ^ and the discussion in the previous section 
we clarified that, in the presence of a nucleus, final-state 
rescattering of the struck small x parton with its rem- 
nants exhausts the similarity between hadronic collisions, 
such as p+A and A+A, and PIS [Icj. In I + A (DIS), 
the multiple interactions of the incoming leptons are sup- 
pressed by powers of a em /a s relative to the struck parton 
scattering. In contrast, in p+A and A+A the initial- and 
final-state scattering of the incoming quarks and gluons 
are equally strong. Nuclear modification, in particular 
jet energy loss associated with the suppression of parti- 
cle production, cannot be neglected. In a model of Su- 
dakov suppression, Ref. 22], it was shown that the effect 
is amplified at forward Feynman x, Xf — * 1- 



A. Evidence for energy loss in cold nuclear matter 

To illustrate the importance of energy loss, we study 
the nuclear modification of hadron production over a 
large range of center of mass energies and momentum 
fractions Xb- In one extreme, at y/s — 200 GeV d+Au 
collisions at RHIC, the STAR collaboration observed a 
factor of 3 suppression of 7r° production relative to the 
binary collision scaled p+p result |l8|. Figure IT^l shows 
that dynamical high-twist shadowing, constrained by the 
DIS data down to values of Bjorken xb ~ 10~ 4 , under- 
predicts the nuclear attenuation by a factor of 2. 

Next, we implement the presently known incoherent 
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FIG. 12: Suppression of 7r° production at y = 4 for d+Au 
collisions at ^/IWn = 200 GeV at RHIC d8 : |. Calculations of 
dynamical shadowing with and without cold nuclear matter 
energy loss are shown. 



limit of energy loss for on-shell partons first derived 
in [l9|. The double differential gluon intensity spectrum 
per scattering reads: 



4 1] 



q? 



du<i 2 k 



k 2 (k- ql ) : 



(22) 



where k is the transverse momentum of the radiative 
gluon, qi is the transverse momentum transfer from the 
medium and 10 is the gluon energy. The two characteris- 
tics of this radiative energy loss that we use are AE oc E 
and AE oc L/X in the case of multiple interactions. Jet 
quenching before the partonic hard scatter implies that in 
p+A and A+A reactions quarks and gluons carry larger 
fraction x a of the nucleon momentum than in p+p colli- 
sions @: 



E A 

<f>q,g/N(x,fif) -> </>q,g/N 



(23) 
(24) 



Fitting Eq. (|23|) to the data, we find k — 0.0175 for mini- 
mum bias reactions implying that an average parton loses 
~ 10% of its energy in a large nucleus such as Au or Pb. 
The variation of k with centrality in the Bertsch-Gunion 
model is given by the nuclear thickness, Ta(6), in a man- 
ner similar to the variation of the nuclear-enhanced power 
corrections. Figure. 1121 shows that such energy loss im- 
plementation gives good quantitative agreement with the 
forward rapidity suppression measured by STAR. 

At the other extreme, the CERN NA35 fixed target 
experiment with y cm — 3 measured hadron production 
in d+Au reactions at y/s^M — 19.4 GeV .'!.'( . Although 
there is a factor ^100 difference in the values of xt,, the 
STAR data were taken at much more forward y than the 
NA35 data, the same rapidity asymmetry is observed as 
at RHIC. We take R y , the ratio of hadron spectra in 
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nuclear modification at three different 
19.4 GeV d+Au collisions, 
theoretical calculations 



J. Right panel 

of charged pion production in p+p and d+Au collisions with 
high-twist shadowing (HTS) and energy loss compared to the 
experiment. 



different rapidity bins defined as 

da h {y) lda h {y ha 







dycPpr I dyd 2 pr 



(25) 



as a measure of cold nuclear matter effects. The left panel 
of Fig. ll3l shows this forward rapidity suppression of neg- 
ative hadrons relative to the baseline y — y cm ~ — 1.5 pro- 
duction cross section. The right panel of Fig. El shows 
the deviation of the data in pt = 0.8 — 1.5 GeV interval 
from the symmetric p+p distribution around y cm . Dy- 
namical shadowing calculations give < 5% effect at 
this energy in this rapidity range. Leading-twist shadow- 
ing parameterizations are also not compatible with the 
asymmetry in the data |l5j . Conversely, the implemen- 
tation of energy loss in cold nuclear matter leads to much 
larger suppression at forward rapidity and significant im- 
provement in the theoretical description of the data in 
Fig. El We conclude that such low energy p+A mea- 
surements [3^ | remove the ambiguity of what nuclear ef- 
fects dominate the forward rapidity suppression in heavy 
ion collisions. 



B. Application of energy loss to D meson 
production in proton-nucleus collisions 



the QCD theory incorporating cold nuclear matter ef- 
fects and the PHENIX measurement of muons coming 
from the decay of light hadrons is achieved. We find 
that the magnitude of the D meson suppression is simi- 
lar to that of pions. It is also comparable with the first 
forward rapidity results on heavy quark nuclear modifi- 
cation in d+Au reactions at RHIC j3(|. In Fig. 1151 we 
also show the result of our calculations for the 7r° — 7r° 
and (D° + D + ) — (D° + D~) mesons pair production at 
large angles. The same rapidity gaps, Ay =y\ — yi, and 
centralities as in Fig. El were used. Nuclear modifica- 
tion of inclusive two particle production is seen to follow 
closely that of single inclusive hadrons. 

There are similarities between the calculations that 
only include resummcd nuclear-enhanced power correc- 
tions, Figs. El and El and the calculations that do 
not neglect energy loss in cold nuclear matter, Figs. ITU 
and 1 151 Both effects are generated through multiple par- 
ton scattering in the medium and lead to suppression of 
the rate of hard scattering. In both cases single inclu- 
sive particle production and large angle di-hadron corre- 
lations are similarly attenuated. Like all nuclear many 
body effects these increase with the centrality of the col- 
lision. The difference in the resulting nuclear modifica- 
tion is that high-twist shadowing arises from the coher- 
ent final-state scattering of the struck small Xb parton 
of the nucleus and disappears as a function of the trans- 
verse momentum. The energy loss considered here arises 
from the initial-state inelastic scattering of the incom- 
ing large x a parton from the proton(deuteron) and leads 
to a suppression which is much more px independent. 
By comparing our results for inclusive one- and inclusive 
two-particle production we see that 



da. 



1 /dyidy 2 dp Tl dp T2 



PT 2 



RdUpTr) 



( N dAu) d,jh NN 2 1 'dyidy 2 dp Tl dpT 2 
da dAu/ d V ld P T l 
( N dAu) ^ nn / d V^ d Pl\ 



Deviation from unity may arise due to specific choices of 
PTi ! Pt 2 and improved theoretical calculations of the par- 
ton color charge, system size and jet energy dependence 
of AE in Eq. (|23[) . Both calculations presented in this 
paper, are consistent with the lack of significant modifi- 
cation of the per trigger yield in d+Au reactions relative 
to p+p [13. 



Having investigated the largest dynamic range of mea- 
surements in proton(deuteron)-nucleus reactions accessi- 
ble to perturbative QCD calculations, we return to D 
meson production and correlations at RHIC. We use the 
same fractional energy loss e = AE/E as in Figs. El 
and El The forward rapidity boost E — tot cosh y\ 3> 
m c allows one to neglect the effect of the charm quark 
mass on AE. Single inclusive 7r° and D° + D + suppres- 
sion at rapidities y\ — 1.25, 2.5 and minimum bias, cen- 
tral and peripheral d+Au collisions at RHIC are shown 
in Fig. El I n this case very good agreement between 



VI. CONCLUSIONS 

In this paper we have systematically studied inclusive 
D meson production and D meson-triggered di-hadron 
yields in p+p and p(d)+A reactions. 

At collider energies, we have found that while DD 
pair production is dominated by heavy quark flavor cre- 
ation processes, such as gluon-gluon fusion gg — > cc and 
light quark-antiquark annihilation qq — > cc, single inclu- 
sive D (and D) meson production is controlled by ra- 
diatively generated charm quarks in the nucleon wave 
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FIG. 14: Nuclear modification of single inclusive 7r° (left 
panels) and D°+D + mesons (right panels) in d+Au collisions 
at RHIC. The same rapidities and centralities as in Fig. llOl are 
shown. Energy loss is taken into account in the theoretical 
calculations. 
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FIG. 15: Nuclear modification of di-hadron correla- 

tions (left panels) and (D° + D + ) — (D° + D~) mesons corre- 
lations (right panels) in d+Au collisions at RHIC. The same 
rapidities and centralities as in Fig. 1111 are shown. Energy 
loss is taken into account in the theoretical calculations. 



function |25l |26j scattering on light quarks and gluons 
eg — ■> eg, cq(q) — > cq{q). Our results for the mecha- 
nism of D meson production are readily testable at RHIC 
and LHC via charm quark triggered two particle mea- 
surements kill Lll . It is a robust prediction of this ap- 
proach that at pt 2 < PTi the D meson triggered away- 
side jet is dominated by hadrons from the fragmentation 
of light quarks and gluons. Only at transverse momenta 
Pt 2 ^ PT t is the associated yield found to be dominated 
by D mesons. Such experimental measurements will be 
able to shed light on the partonic processes responsible 
for D meson production and also constrain the hardness 
of the non-perturbative heavy quark fragmentation func- 
tions p^. Our results are also suggestive that inclusive 
heavy meson cross sections may be larger than the ones 
obtained in perturbative expansions that always treat c 
and b quarks as "heavy" |24j |. 

In p(d)+A reactions we calculated and resummed the 
nuclear-enhanced coherent power corrections 0, 0] 
to single inclusive D meson production and D meson- 
triggered di-hadron yields. Improving upon our previous 
work [h}, we established here the effect of the heavy 
quark mass on the power corrections and found it to be 
relatively small for charm quarks. We also extended the 
study of coherent multiple scattering to all s-, t- and it- 
channels. While forward rapidity at RHIC is dominated 
by DIS-like i-channel final-state interactions other kine- 
matic domains and other particle production processes, 
such as the Drell-Yan, are controlled by different chan- 



nels. In this paper, by direct calculation, we demon- 
strated that initial-state s-channel multiple interactions 
lead to an enhancement in the cross sections in clear con- 
trast to final-state suppression. We have thus further cor- 
roborated the conclusion that enhancement/suppression 
effects in reactions with nuclei are dynamically generated 
in the collision d H [H |U E3 We find that they are 
not universal (are process dependent) and cannot be sep- 
arated in the parton distribution or fragmentation func- 
tions. 

Numerically, we found that dynamical high-twist shad- 
owing can suppress single inclusive D meson production 
by as much as it does light hadrons. Similar results were 
obtained for di-hadron production with light and heavy 
quark triggers in the rapidity range y\ = 1 — 4 studied 
at RHIC. The nuclear modification at forward rapidity 
and small pr ~ 1 GeV can be substantial, ~ 30 — 50%. 
At high pt, however, it disappears quickly oc \jp\ and 
cannot fully explain the nuclear suppression measured 
in d+Au reactions by PHENIX [HlEi and STAR [3. 
Our calculations are fully constrained by existing DIS 
data |12| and provide an upper limit on nuclear modifi- 
cation from dynamical high-twist shadowing. Deviations 
between the calculations and the existing data are indica- 
tive of inelastic scattering in nuclei j2^| . 

To determine the additional dynamical effects that lead 
to forward-backward hadron production asymmetry in 
proton-nucleus collisions in a model-independent way, we 
investigated low energy d+Au collisions [33| where coher- 
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ent scattering does not play a role. Suppression of for- 
ward rapidity hadrons, similar to that observed at RHIC, 
is found and is consistent with jet energy loss in cold 
nuclear matter. We have shown that, based on initial- 
state partonic implementation of energy loss @ similar 
to the established incorporation of final-state AE 
the forward rapidity suppression can be described over 
the widest dynamic range in p(d)+A collisions accessi- 
ble to perturbative QCD analysis. Specifically, we im- 
plemented radiative energy loss of the type derived by 
Bertsch and Gunion an d found that if the average 
parton loses ~ 10% of its energy in minimum bias re- 
actions with nuclei such as Au and Pb, this can explain 
the suppression of forward hadron production from SPS 
to RHIC. Inclusion of cold nuclear matter energy loss 
together with coherent power corrections results in sup- 
pression of both single inclusive hadron production and 
hadron triggered inclusive two particle yields (includ- 
ing D mesons and DD pair production) which are fully 
consistent with the forward rapidity hadron attenuation 
measurements at RHIC llll . 

The implication of such energy loss for extracting the 
scale of momentum transfer in cold nuclear matter mer- 
its further investigation. Of particular interest are an- 
alytic solutions for initial-state gluon bremsstrahlung, 
including coherence effects to all orders in the correla- 
tion between multiple scattering centers. In the future, 



These can be studied numerically and implemented in 
the pQCD formalism with a degree of sophistication that 
will match the well-developed final-state energy loss phe- 
nomenology p|. 

In summary, we have presented the first perturba- 
tive QCD calculation on heavy meson triggered di- 
hadron yields and showed how it can provide informa- 
tion on the D meson production mechanism and the non- 
perturbative fragmentation of heavy quarks. We identi- 
fied the nuclear effects that dominate the cold nuclear 
matter attenuation of D meson and light hadron pro- 
duction at forward rapidities, which has drawn consid- 
erable interest in its own right. This work also provides 
the baseline for precision QGP tomography using heavy 
quarks as probes of the plasma away from midrapidity 
where new experimental capabilities at RHIC and the 
LHC are expected to become available. 
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APPENDIX A: FLAVOR CREATION AND FLAVOR EXCITATION DIAGRAMS 



We emphasize that our need for the matrix elements calculated here is our interest in matching the m c = charm 
quarks from the parton distribution functions to the m c ^ in the final state. The internal heavy quark propagators 
are treated as massive. 

The dominant mode of charm-anticharm di-jet production and subsequent fragmentation into DD mesons comes 




FIG. 16: (a) Light quark-antiquark annihilation into a massive cc final state, (b) t— channel scattering contribution for the 
case of non-negligible initial (anti)charm parton distributions. Note that for cc — > cc both diagrams contribute. Time increases 
from left to right in all diagrams and arrows indicate momentum flow. 
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FIG. 17: (a) Gluon fusion via the s— channel into a cc final state, (b) Ghost- ant ighost annihilation that cancels the contribution 
from the unphysical gluon polarization states in the covariant gauge, (c) t— channel and (d) the crossed u— channel charm 
exchange diagram. 



from gluon fusion and light quark- antiquark annihilation. The cc creation diagrams from q = (u, d, s) quarks are 
shown in Fig. I16f a~). An additional t— channel scattering diagram contributes in Fig. I16f b'l in the case of non- 
negligible (anti)charm distributions in the nucleon wave-function. In all figures the arrows indicate momentum flow. 
The matrix element corresponding to Fig. Eta) reads: 



iM 



(a) 

qq^cc 



(pi + P2) 2 + 

For identical initial- and final-state quark flavors, the t— channel must also be considered: 

ib 



iM 



(b) 
cc — >cc 



u(h) (-ifl,yT°)u(pi) 7 . v(p 2 ) (^g sl »T b ) v(k 2 ) 

[pi - ki) z +ie 



(Al) 



(A2) 



The relative "— " between Eqs. (|A1|) and (|A2|) comes from Wick's theorem. The dominant contribution to the cc cross 
section comes from the diagrams with gluonic initial states: 



iM 



(a) 

gg^cc 



-igw 



(pi+p 2 y + ie 

(9sf abc ) [«T (Pi - P2) X ' + g vX ' (2 P 2 + PiT + 5 AV (~2 Pl - p 2 y] e x - ( Pl )e^ (pa) 



(A3) 



iM 



(c) 

gg^cc 



u(h) (-i;, 7 T)i 7 ' (tl (-WZ*) v(k 2 ) e^( Pl )eiHp 2 ) 



(A4) 



hjW - .-/u.u_ s „~mt*\.- 7-(Pi-fe2)+m c ; , / \,,„. w a j 



iM^ = S(fc 1 )(~ l .g s7 T b )z „ ., . 

(pi — k 2 y — m£ + it 



{-ig s ^T a )v{k 2 )e^(p x )e^{p 2 ) 



(A5) 
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FIG. 18: (a) Heavy quark on light quark scattering via t— channel gluon exchange. Diagrams for scattering of anticharm and/or 
light antiquarks are similar. 



In the covariant gauge used in this calculation, there are contributions at the cross section level coming from the 
unphysical gluon polarizations. The diagram M.^ in Fig. lltjf a) is an example of where Fadeev-Popov ghost-antighost 
annihilation diagrams, two of them with fj and rj interchanged, have to be taken into account: 

iM^ cE = u(kx) (-ig s rT b ) v(k 2 ) - ~^ 2 V {g s f abc ){- P i) ■ (A6) 

(P1+P2) +i£ 

We recall that a closed ghost loop comes with a "— " sign and the color factor is the same as in jA^"- 1 ! 2 . 

For the case of single inclusive D meson production, (anti)charm scattering on light quarks and gluons must also 
be included. We illustrate the calculation for quarks, as shown in Fig. Results for antiquarks are similar. Only 
one diagram contributes to heavy on light quark scattering: 

iM^ cq = -u(h) (-ig s ^T a ) u{pi) ~y^ 2 ab + . u(k 2 ) (-ig s YT b ) u( P2 ) . (A7) 

Figure 1191 shows the contribution of charm quark scattering on gluons. For t— channel gluon exchange, the matrix 
element, , is: 



iM%Lcg = S(*i)(-i«7.7 A r c )«(Pi): 



(fci - pi) 2 +ie 

(g s f abc ) [<T (p 2 + hf + g^' (p 2 ~ 2k 2 f + g x '»{k 2 - 2 P2 y] # (p 2 )e^ (k 2 ) , (A8) 
where the unphysical gluon polarizations are again canceled by the ghost contribution: 

iM^ cg = u(h) (~zg sl x T c ) u( Pl ) " tgAA 2 ' . (g s f abc )(- P2 ) x ' . (A9) 
J J (fci — p\) + ie 

The other diagrams that contribute are it— channel and s— channel virtual quark exchange, A4^ and A4^ respectively, 
given for the example of charm by: 

iM^ cg = u(h) Hg s7 ^)» 7 ' (Pl ~ fc2)+ 2 r " c {-ig sl »T b ) u(pr)e^(p 2 )el x *(k 2 ) , (A10) 
J a (Pi — k 2 ) z — mf, + ie p 



iMi d g l cg = fifo) {-ig^T^i 1 ^ 1 t k2) + m j _ c . i-igs^T*) ufrOefr (p^fo) , (All) 

[Ki -f- K 2 ) m c IE 

Next, we recall the definition of the Mandelstam invariants: 

s = ( P i +P2) 2 = (h + k 2 ) 2 , 
i = {Pi ~ h) 2 = (p 2 ~ k 2 f , 

u = (pi - k 2 f = ( P2 - h) 2 . (A12) 
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FIG. 19: (a) Charm quark on gluon scattering via t— channel gluon exchange, (b) Charm-ghost scattering that cancels the 
contribution from the unphysical gluon polarization states in the covariant gauge, (c) Charm quark on gluon scattering via 
virtual heavy quark exchange in the u— channel, (d) The corresponding s— channel diagram. Anticharm on gluon scattering is 
given by a similar set of diagrams. 



In simplifying the result for the squared matrix elements, for charm-anticharm production we have s + t + u = 2m 2 
and for charm on light parton scattering we have s + t + u — m 2 . We have taken all massless initial-state partons. 
The spin (polarization) and color averages (• • • ) are defined in the usual way: 



(...) 



spin 



olor 



- E 

s=±l/2 

- E 



a=l---N a 



quarks , 
— quarks 



' ) polar 



ization 



A=1.2 



(N- 



gluons 



— E 



a=l— JV?-l 



gluons , 



(A13) 
(A14) 



where N c — 3 is the number of colors. Summation over the final-state spin (polarization) and color is understood. 

In calculating the color factors we used factorization of the color algebra and the Lorentz structure in the Feynman 
graphs. A diagrammatic approach in non-Abelian gauge theories is employed in the calculation below. This is the 
same approach that was used in the evaluation of the significantly more involved color factors [23, 0] associated with 
medium-induced gluon bremsstrahlung. The calculation of the squared amplitudes, Eqs. (|A1(I - (|A10(I is straightforward 
though somewhat tedious in the presence of non-zero parton mass for final-state heavy quarks, XL=±i/2 u s (p)u s (p) = 
7 -p + m, J2 s =±i/2 v s (p)v s (p) = 7 -p — m. In the squared matrix elements given below, the color factors are separated 
in angle ( ) brackets and a "— " sign inside denotes a destructive color interference. 

For light quark annihilation into a massive cc final state: 



(\M q ^ c5 \ 2 ) 



Amis 



2mt 



(A15) 



For charm quark- antiquark scattering by t— channel gluon exchange and s— channel annihilation into a massive cc pair 
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the result is: 



2\ 2. 9 



2\ 2ft 



9/ i 2 



(A16) 



Finally, the result for gluon fusion into a heavy charm-anticharm pair, which includes the corresponding t— channel 
and u— channel scattering, is: 



12/ (t-m 2 ) 2 
3 ^ ^ 4 



1 \ 2 3 . 



32/ s(i-m 2 ) 



(f 2 + m 2 s-2m^ + m^) + (- — 



-— (tu + m 2 s — 2m 2 u — imt) 

(U - TO 2 ) 2 V y 

( 2 2 2 4 \ 

it + m c s — 2m c u + m c ) 



32 / s(u — m 2 ) 



96 / (t — m 2 ) (u — m 2 ) 



(m 2 s — 4m*) 



(A17) 



We turn next to the additional contributions to single inclusive charm production. The scattering of heavy on light 
(anti)quarks via t— channel gluon exchange yields: 



\\M,„ -,„| 2 } = ( 2 -y-^(s 2 + u 2 + m li- m i) 



(A18) 



For the case of charm quarks scattering on gluons, we obtain: 
(\M cg ^ cg \ 2 ) = +(^^(Ai 2 ~Asu-m 2 c .s-3m 2 J-mi) 



9 / (u — m 2 ) 2 
1\ 2ftf 



{su + 2m 2 u — m 2 s) — / — 



2 5s 4 



9 / (s — m 2 )- 



4/ t{u — m 2 ) 



36/ (s — m 2 )(u — m 2 ) 



(2u 2 5m 2 c u + mli-mi) + (l) 

\4/ <(s — rnjfj 



2s — 5m r s + m r t — m,,\ 



(A19) 



We can partially check all results by taking the massless, m c — > 0, limit. If we denote by (q)q the light u,d,s 
(anti)quarks, we recover the familiar form of the squared matrix elements: 



ro lim(|A^ c5 | 2 ) = (K„ y f) = lf(^- 



4ft* ft 2 s 2 +u 2 2u z 



lirn <|A4 C ^ CC -| 2 ) = (\M q ^ qq \ 2 ) = 



f 2 



3 si 



Jim <|M c ^ ct( | 2 ) = <|M ? ^| 2 > = ^(^ 



(A20) 
(A21) 
(A22) 
(A23) 



lim (IMcj^cg 



l 2 > = (IX»-»I J > = si 



9 \ u s 



For comparison, see Ref. @|. 



(A24) 
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APPENDIX B: FACTORIZATION APPROACH TO COHERENT MULTIPLE SCATTERING 



In proton-nucleus reactions at RHIC the negative lightcone direction is assigned to the nucleus. Thus, the proton 
has momentum pa — [p\~, 0, OjJ and the nucleus has momentum Pb = [0, P B , 0j_]. We again remind the reader that 
11 B" and "6" are just labels for the target nucleus with mass number A. For example, pb ~ PbI A is the momentum 
per nucleon. Here, we neglect the Fermi momentum and the binding energy of the nucleons in the nucleus. For 
momentum exchange q^ = (p c /zi — XaPA)^ we rotate to a frame where parton "d" interacts head on with the 
remnants of the nucleus. Let = [1, 0, OjJ, n M = [0, 1, 0j_] specify the "+" and "— " lightcone directions, respectively. 
In the equivalent Breit frame 



-x b p B i 



2x b p B 



n , 



t . 



(Bl) 



Here, x b can be derived from the requirement that when the scattered parton "d" of momentum xp B n^ 
shell. The t — oo cut constrains the momentum fraction from the nucleus to the known solution, Eq. ©, 



q^ is on 



S{(q + xpsf - m|) cx 5(x ~ x b ) 



x b 



l + m3/(-t) 



(B2) 



We note that above, u d" is the label, not the flavor, of the outgoing parton. In particular when "d" is a c quark we 
will have = m c = 1.3 GeV. 

Under the standard convention about the "+" and "— " lightcone directions at RHIC, it is convenient to use the 
A ■ fi = A~ = gauge. Such a choice does not affect any of the physical quantities that appear in the final result. 
For example, the parton distributions, the scale of higher twist and the hard scattering partonic cross sections and 
the Mandelstam variables are all individually boost and gauge invariant. We consider the final-state interactions of 
parton "cf with the gluons from the nucleus with A^{y^) w A ± {y^) shown in Fig. [201 Using 



dA^(y+) dA-(y+) 



9Vi 



lg [A^(y+),A-(y+)] 



dvf 



(B3) 



the gluon fields can be represented as components of the field strength tensor. The initial-state parton carries 
momentum fraction x$ and 2n gluons carry momentum fractions Xi — a^-i, at positions yf, and fractions Xi — Xi, at 
positions yf, of the momenta pb = Pb/A of their corresponding nucleons. Such choice of the momentum transfers 
ensures the flow indicated in Fig. 1201 and a simple representation of the momentum fraction in each of the propagators 




X S I N X > 



A 17 A 17 



A lf A 1 ' 




t = 00 



A 



FIG. 20: Multiple coherent scattering of the outgoing partons in proton-nucleus reactions in the t-channel. We have denoted 
by and -j — > the long distance and contact propagators, respectively. We have indicated the lightcone positions, e.g. yf. 
Arcs show the momentum routing. Gluons are dominated by their transverse components and the fields, e.g. A ± '(yf), are 
shown before the conversion to the field strength tensor and before the factorization of the color and Lorentz indexes. 
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along the scattered parton line. The part of the diagram associated with the non-perturbative matrix element reads: 



I I 2ir 2ir 

x (p B \o init A ± (y+)A^(y+) ■ ■ ■ A ± (yf)A^ (yf) ■ ■ ■ A ± (y+)A^ (y+)\p B ) 



2tt 2ir i(xi — Xi—i — ie) i(xi — Xi — ie) 
„ , pg i oin u F -± {y + )Fli ~+y.. F -± {y + )Fl ^ Pb) _ (B4) 



In Eq. I|B4|) we have carried out an integration by parts assuming that the gluon fields vanish sufficiently fast at 
— > oo, — > oo. The choice of — ie is convenient for guiding the evaluation of the remaining momentum 
fraction integrals [l2T ] but can be verified by carrying out twist-4 calculations in the lightcone and covariant gauges 
and comparing the final result. Depending on the type of initial-state parton, the operator 



0^=^(0)2-^ 



1 



=F- x (0)FZ(y+) . 



(B5) 



To factorize the twist-2 + 2n non-perturbative matrix element we need to decouple the color and Lorentz indexes, 
which are not shown explicitly in the equation above. We use Fiertz transformations and identify the numerical 
factors associated with the diagonal terms that do not change the quantum number, for example color, of the parton 
after a pairwise gluon exchange. We obtain —1/2 for the Lorentz part and l/2N c — 1/6 or N C /(N^ — 1) = 3/8 
for a final-state interacting quark or gluon, respectively. After this factorization of indexes, we obtain pairwise 
contractions F~- Ll (y~f)F^ Li (y~f) in the non-perturbative matrix element, where average over color is implicit, and 
— h{(j±)ijPI • ■ ■7" Ll in the hard scattering matrix element, which we can evaluate perturbatively in QCD. 

Let us now focus on the propagator structure along the outgoing parton "d" . For definiteness we consider a heavy 
quark. With net momentum transfer Xip^ from the nucleus, the quark propagator can be represented as follows 



T>(xi,Xb,—i, m d ) 



Xb 



PbT 



Xb \ (7 " P + m d) 



Xb 



(B6) 



In Eq. i|B6|l the Xi-independent vector p is on shell 



(xb - x b )p B n' J - + 



2x b p B 



(B7) 



We emphasize that the propagators along the final-state parton line have two distinct parts |12t Il3j . The first one 
does not have a pole in terms of the yet undetermined momentum XiPg. In the Fourier space of the conjugate variable 
y + , the interactions that such a propagator separates can be collapsed to a single point, oc S(Ay + ). Therefore, the 
short-distance two-gluon exchange can be evaluated within the same nucleon state. On the other hand, the second 
part of the propagator has a pole at Xi = x b , which puts the propagating parton on shell and ensures a long lifetime. 
In Fourier space, the subsequent interactions separated by this propagator cx 6(Ay + ) can be long distance and are 
only limited by the nuclear size ~ roA 1 / 3 , where A is the nuclear target mass number. In our notation the t — 00 
on-shell cut line through parton "d" reads 



C(xi,Xb,-i, m d ) 



2n [ -| ) (p • 7 + m d )S(xi - x b ) 



(B8) 



Naively, the sequence of propagators separated by the j 1 - matrices may yield an exponentially large number of terms 
per single digram, similar to the one in Fig. 1201 However, using the equations of motion for a massive particle, the 
commutation relations for the Dirac algebra and (7 + ) 2 = (7~) 2 = 0, we can simplify a typical sequence of propagators 
along scattered parton "e?" line as follows: 



(p ■ 7 + m d ) 
(p ■ 7 + m d ) 



9t 
2N C 

2ira s 



Xb \ 7 T 7~ 



Xb 

-t 



1 



x b + ie 



j-p + m d , ± 

~ — - 7 Pb7 

Xi — Xb + ie ' ( 



7 ■ p + m d 
Xi-i - Xb + '1 



(B9) 
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where the boundary condition is set by the cut, Eq. I|B8|) . and we have taken into account that the subsequent 
scattering has similar propagator structure. This is the first significant simplification, which makes it possible to 
calculate and resum the multiple final-state interactions of the struck parton. 

The second important simplification comes from a model [T3| for the high-twist matrix elements in Eq. I|B4|) . The 
natural long and short distance separation of soft gluon interactions along the trajectory of the struck parton justifies 
the decomposition of the multi-local, multi-parton, non-perturbative matrix element in Eq. (|B4[1 in a diagonal basis 
of individual nucleon states. We focus on the QCD aspects of the calculation and not the ones related to the nuclear 
structure. To illustrate our idea, we adopt a simplified model of a nucleus of constant lab frame number density: 

P(r) = £ = 4^ > ( B1 °) 

where r = 1.2 fm is the nucleon radius. Each nucleon carries an equal fraction of the momentum of the nucleus, ps ~ 
Pb/A, an approximation which is valid so long as the energy per nucleon is significantly larger than the binding energy, 
e ~ 8 MeV, and the momentum per nucleon us significantly larger than the Fermi momentum, kp w 250 MeV. If this 
condition is violated, the off-diagonal elements in a density matrix that represent the correlation between different 
nucleons might contribute significantly. With the normalization of the momentum states (p' \ p) — 2E p (2ir) 3 8 3 (p — p'), 
the decomposition of the matrix element reads: 



O mlt Y[F-^(y+)FZ(yt) Pb) « A {pbP^pb) \{ ^(^^(^1(^)1^} 

3 



2E P 



= A{p B \O init \p B )]l— 3 (p B \F-^(yt)F-(yt)\ P B) ■ (Bll) 

From Eq. i|Bll() it follows that, in the absence of final-state coherent interactions, i.e. neglecting the Y\a—i two- 
gluon correlations, the minimum bias, inclusive, deeply inelastic scattering structure functions, and the minimum bias 
differential hadron production cross section in proton-nucleus reactions, both scale with the number of the nucleons, 
A, in the nucleus. The expectation value of the bifocal quark or gluon operator O m%t enters the definition of the 
leading-twist quark parton distribution function in the nucleon and is always present as a first term in the twist 
expansion series. 

Next, we carry out the integrations Oi==i J dxidxi. To determine the scale of high-twist corrections per nucleon, 
we recall that the integration over the single pole x l in Eq. (|B4|I yields terms proportional to 2ir9(yf > yf) and 
2n9(yf < yf)- If the two gluon exchange is to the left of the t = oo cut we keep the term cx 0(yf > yf) and if 
the two gluon exchange is to the right of the t — oo cut, we keep the term oc 9(yf < yf). This ensures that time 
grows toward t = oo and that we keep the dominant nuclear-size-enhanced contribution. Inserting 1 = exp(i0p £ >y + ) 
we identify 

(2tt) / ^e*°^ + — (p B \FZ{yt)F-\yt)\PB)6{yt > or <yf) = (2tt) lim ^xG(x) (B12) 
J Zn p B x-^a i 

by direct comparison with the definition of the parton distribution functions [T^ . Note that the lightcone positions 
yf 1 yf can always be shifted due to translational invariance. The ^-function yields a factor of 1/2 which can be 
thought of alternatively as arising from the symmetry in the legs of the two-gluon ladder at the interaction point. We 
emphasize that the scale of high-twist corrections depends on the soft x — > limit of the gluon momentum fraction 
in the nucleon. While finite resolution effects may alter this result to some small but finite value, this is the value 
of neither the Bjorken x nor the momentum fraction fixed by the hard scattering in hadronic reactions, Eq. ©. As 
emphasized below Eq. (|2J), we explicitly keep track of the A n l 3 enhancement of the multiple soft scattering. We will 
show below that taking all possible cuts in a diagram with n 2-gluon ladders gives contribution oc (l/nl)A n / 3 . We 
first identify the prefactors associated with such an exponentiating series. A dimensionless equivalent of an average 
path length of a struck parton through the nucleus is given by 

/ / dy+p B 6(yf)) = -m N r A l l* , (B13) 

\J I Volume 4 

where the (• ■ ■ ) volume is carried out for uniform density, Eq. I|B10|) . In a sequence of interactions where the two- 
gluon ladder contribution does not vanish, the 7 + 7~/2 matrices give a numerical factor of unity; for example, 
Tr 7 -( 7+7 -/ 2 )„ 7+ . . . = Xr 7 - 7 + . . . < Collecting all factor from Eqs. (BUJ, (|BTT|l and ETgL we obtain the pa- 
rameter that controls the strength of the high- twist corrections on a single nucleon (A=l) [lOL Ha. Il3| 
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FIG. 21: Reduced Feynman diagram to the left and right of the t = oo cut. The Feynman rules for the effective scalar 
interaction are given in Eqs. 1B15II and (113161 . 



As we mentioned above, this does not include the A 1 / 3 nuclear size enhancement, nor the derived factor of 1/nl which 
can be intuitively understood from path ordered interactions pol | . 

It is useful to define effective Feynman rules, corresponding to the diagrams shown in Fig. 1211 for the final-state 
soft scattering of the struck parton: 



a i{xi— Xi_i)Ai 



27T Xi — Xi-i — ie Xi-\ — Xf, + ie 
dX, e «(^-^-i)A< -i 



2ir Xi — Xi-i — ie X{ 



x b 



ic 



left of t = oo , 



right of t = oo 



(B15) 
(B16) 



In Eqs. I|B15|) and l|B16j) we introduced the notation \i — Pgy^ ■ The difference from our result for DIS, given in 
Ref. [12j, is that we have already performed the density matrix decomposition of the high-twist matrix element. Here, 
fb replaces Xb and depends on the heavy quark mass, and Q 2 — ► — t. Numerical factors, including the contribution 
to the non- vanishing traces from the 7-matrices, are collected in the scale of higher twist, £ 2 , Eq. l|B14jl . 

The last part remaining is to demonstrate that the dimensionless integrals can be performed yielding the claimed 
A 1 / 3 amplification of the higher-twist contribution through the exponentiation of the power series. The twist-2 + 2n 
term reads: 



J2 [ dx d Xl ---dx n [ ^^...^ e ^°±tp B \0™\p B 
J 2?r 2tt 2?r p B 



n 



3=1 J 



Xb + ie Xj — Xj-\ — ie 



— 8{xi - x b ) 



n 

k=i+l 



i(x fc -x fe _ 1 )A fc 



Xk - x b - ie x k - x fe _i - ie 



(B17) 



In Eq. (|B17|I . the appearance of is associated with the conversion of to Ao- There may also be an additional 
function, H(xq), which contains smooth xq dependence. The n = first term corresponds to there being no final-state 
interactions. 

The integrals to the left and right of the t = 00 cut, which sets taken independently. The integrals over 

the Xj, Xk poles can be carried out by closing the contour in the upper or lower half plane. We note that we have 
multiple poles. These produce 9- functions which lead path ordering of the d\ a integration and powers of X a — Xp. 
Carrying out the J dXi ■ ■ ■ dX n , we find 



dXy 
2ir ' 



EbXo — (p B \O mlt \p B ) V 
P B ^ 



i=0 



;(A-A ) ; 



(n-i)\ 



:(A) 5 



/ ^e ix > x ° — ( P b\0" m \pb) i (*AA )" = / ^e^vt (p B \O init \p 

J 2tt p B n! n! dx b l [ J 2tt 



(B18) 



In Eq. llBTSll . A = A 1 / 3 A and we have dropped a term 0(Xq/2) which is suppressed by the large A 1 / 3 nuclear size 
relative to AAq. In the presence of an additional function of xq, for example H(xq) — \A4 a b-+cd\ 2 /%o, the multiple 
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poles to the left of the t — oo cut lead to a series of derivatives on H(xq). The result for the modified Eq. I|B17() . 
where we have also included the (istbS, 2 / (—i)) n power dependence, reads: 

£ (^§y 1 / dx Q dx x ...dx n j ^^...^ e ^°-L(p B \G^\p B ) H(x ) 



n 



— — Xj-i — x\y + i€ Xj — Xj^i — ie 



8(xi - x b ) Y[ 



i(a; fc -3: fc _ 1 )A fc 



k=i+l 



Xk- Xb - it x k - - ^e 



n\ \ Xh -t ) dxl 



H(x b ) 



2ir 



e^oPaVo (p B \O mlt \p B } 



(B19) 



Note that all infrared singularities, naively present on the hard part, give a finite contribution after taking all possible 
cuts in the diagram, Fig. [3U1 

In studying the modification to the differential cross sections due to nuclear-enhanced power corrections, we can 
either take the scale defined in Eq. I|B14(I amplified by the nuclear size oc A 1 / 3 and a similar effect for the proton with 
A = 1, or assume that, for the case of the "elementary" nucleon, high- twist effects have been taken into account in 
the PDFs and additional interactions in the nucleus are oc A 1 / 3 — 1. The extraction of the scale of power corrections, 
£ 2 = 0.09 — 0.12 GeV 2 , for quark scattering, from comparison to the shadowing in the nuclear structure function 
F^(x,Q 2 ) was performed using the second assumption. Improved analysis may yield somewhat different values 
for £ 2 . Finally, we recall that in the case of final-state gluon scattering, the color singlet coupling, N c /(Nq — 1) = 3/8, 
amplifies the scale of higher-twist effects by a factor of Ca/Cf = 9/4 [To| . 

To complete the summation of all-twist, nuclear-enhanced power corrections, we identify the additional function 
H(xb) — | M a b-^cd\ 2 / Xb from Eqs. 10 and i|12|l . In Eq. ljB17J) . it multiplies the leading-twist parton distribution 
(t>b/N( x b)- Thus, if we define fib] F a b^ c d{x b ) = <pb/N( x b)\Mab~*cd\ 2 /xb- Implementing the normalization to no nuclear 
effect on a single nucleon as follows, A 1 / 3 — > A 1 / 3 — 1, and having allowed for both massless and massive partons 
(quarks or gluons) we find: 



v 2 * 1 ( £ 2 U 1/3 -l)\" d n 
F ab ^ c d{x b ) — I I Xb j ] -j-^F a b-, c d{Xb) = exp 



dxl 



£2(^1/3 _ 

F a b~*cd [ Xb + XbCd ? I = F a b^cd [ Xb 



£2(^1/3 _!) d 
Xb — 



l + C d 



—t dxb 



F a b^cd(xb) 



-t + ml 



(B20) 



In Eq. (|B20(I Ca = 1 (9/4) for quarks (gluons), respectively. When md -> we recover the result of ^{j- The mass 
here reduces the effect of power corrections and compensates for the fact that Xb already has the physical mass of 
parton "d" factored in. 
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